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Abstract 

The fusion hindrance in heavy ion colhsions is studied in the framework of the two-center hquid 
drop model. It appears that the neck and the radial degrees of freedom might both be hampered 
by an inner potential barrier on their path between the contact configuration to the compound 
nucleus. Heavy ion reactions with and without the two kinds of fusion hindrance are classified 
through systematic calculations. It is found that the number of reactions without radial fusion 
hindrance is much smaller than that without neck fusion hindrance, and for both kinds of fusion 
hindrance the number of reactions without fusion hindrance at small mass-asymmetry parameter 
a is smaller than that at large a. In the formation of a given compound nucleus, if a reaction 
with etc is not hindered, then other reactions with a > Uc are also not hindered as it is well known 
experimentally. 

PACS numbers: 25.70.Gh, 25.70.Jj, 27.90. +b. 
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I. INTRODUCTION 



The fusion hindrance that appears in heavy ion reactions has been known for many years. 
Now the mechanism of fusion hindrance is gradually understood as to be due to an extra in- 
ternal barrier between the touching configuration and the compound shape after overcoming 
the Coulomb barrier [l-3|. However what causes the internal barrier is a difficult problem. 
One of the hypotheses is that the internal barrier could be thought as the conditional saddle 



d be attributed to an effective barrier due 



point in the liquid-drop potential as well as cou 
to the dissipation of the incident kinetic energy 

Since there are two barriers for the fusion, its theoretical description is divided into two 
consecutive steps: one is from infinity to the contact configuration passing over the Coulomb 
barrier, and another one is from the contact point to the compound state overcoming the 
internal barrier. This so-called two-step model explains the ext ra-p ush energy and further- 



more gives an energy dependent fusion cross sections , [lO|. The measured residue 

cross section can be calculated by plugging a statistical evaporation model on the two-step 
model. 

To evaluate the amplitude of the hindrance to the fusion, the size of the inner barrier is 
not enough. The dissipation mechanism also plays an important role and one has to solve 
dynamically the diffusion over this barrier. This was done analytically in one dimensional 



case assuming a parabolic barrier ll|, ll2|. Here, we will focus our attention on the appear- 



ance of the hindrance to the fusion that is due to the saddle point height measured from 
the energy of the touching configuration formed by the projectile and target nuclei of the 
incident channel. In other words, there is no hindrance for cases without internal barrier. 
Therefore, the border between the normal and the hindered fusion is given by the condition 
that the contact configuration be on the top of the saddle point. 

For mass-symmetric reactions, the fusion path between the touching configuration and 
the compound shape can be easily determined since the system evolves only along the radial 
degree of freedom. In order to determine if the fusion hindrance occurred for a given reaction, 
we simply compare the position of the barrier and the touching configuration: the fusion 
exists while the barrier is located inner side of the touching point and vice versa. Ref. fisj ] 
systematically analyzed different symmetric reactions and drew a clear border between the 
two kinds of reactions (with and without fusion hindrance). 
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How is the situation for mass-asymmetric cases? Generally this problem becomes difficult 
because the system crosses a ridgeline, not a saddle. It is then difficult to determine an 
effective reaction coordinate. However there is a way to avoid the difficulty: let the reactions 
evolve starting from the contact point on the potential surface without any random force, 
and then check if the compound state is reached. If the composite system evolves toward the 
compound direction, it means that there is no barrier between contact point and compound 
state, and thus the fusion hindrance does not exist. Otherwise, the fusion hindrance occurs. 
With this method, the current paper studies the occurrence of the fusion hindrance for 
asymmetric reactions. 
Refs. 
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18| show that the neck of the amalgamated system also plays an important role 
in the hindrance to the fusion. Thus, one of the related questions is about the evolution of 
neck degree of freedom: Is there also a fusion hindrance in the neck evolution, like for radial 
evolution? 

We will answer to these questions in the present paper. 
II. PARAMETERIZATION OF THE COMPOSITE SYSTEM 



;he amalgamated system. In this 
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15| based on three important 



There are several ways to parameterize the shape of 
paper we use the so-called two-center parameterization 
parameters which are: the distance between two centers z, the mass-asymmetry parameter 
a, and the neck parameter e. The dimensionless parameter z is defined as follows: 

z = R/Ro, (1) 

where R denotes the distance between two centers of the harmonic potentials, and Rq the 
radius of the spherical compound nucleus. The mass-asymmetry parameter is defined as 
usual, 

where Ai and A2 are the mass numbers of the colliding nuclei. The neck parameter e is 
defined by the ratio of the smoothed height at the connection point of the two harmonic 
potentials and that of spiked potential. In this description, nuclear shape is defined by 
equipotential surfaces with a constant volume. For example, e = 1 means no correction, 
i.e., complete di-nucleus shape, while e = means no spike, i.e., flatly connected potential. 



which describes highly deformed mono-nucleus. Thus, the neck describes shape evolution of 
the compound system from di-nucleus to mono-nucleus. The initial value of 2;, a and e, for 
a given reaction at the touching point is 

^1/3 ^ ^1/3 

"° = (X + A)V3' (3) 

= 1, (5) 

respectively. Here (At) stands for the mass number of the projectile (target). With this 
parameterization, the liquid drop model (LDM) potential is adopted to calculate the energy 
of the touching di-nucleus. 

These three parameters are connected by the potential, the inertia and friction tensors. 
For the sake of comprehension, we will study the neck and the radial degrees of freedom 
separately. 



III. FUSION HINDRANCE IN THE NECK EVOLUTION 

In this section we focus our attention on the potential landscape seen by the neck pa- 
rameter. Let us recall that in the two-step model, the fusion hindrance is described as the 
additional barrier hampering the evolution from the contact point to the compound nucleus. 

In order to show the direction of evolution of the neck degree of freedom, the LDM po- 
tential is plotted as a function of e in Fig. [H where the different lines correspond to different 
distances between the projectile and target. The figure shows two different directions of 
evolution: (i) at small z i^z = 1.3 or z = L5) the potential decreases with reducing neck 
parameter (thicker neck) which means that at contact the neck becomes thick by the driv- 
ing effect of the potential, (ii) at large z {^z = 1.8 or z = 1.9), the potential increases with 
reducing neck parameter around e ~ 1 and a bump exists at a critical value of e. This 
indicates that an extra energy should be provided to reach the thick neck, i.e. to form a 
mo no- nucleus. In other words, the system is hindered to form a compound nucleus. In fact, 
the frontier between the two cases can be determined by the slope of the potential at the 
point e = \: 

„ dV 
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FIG. 1: LDM potential as a function of the neck parameter e for ^^Kr+^*^^Pb at various relative 
distances. 
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FIG. 2: LDM potential for ^^°Pd + ^^°Pd. (a) LDM potential as a function of e; (b) Components 
of the LDM potential, surface potential (above zero) and Coulomb potential (below zero), as a 
function of e; (c) Di-nuclear shapes at e = 1. The dashed, solid and dash-dotted lines are for 
z = 1.4, 1.67, 1.8, respectively. The slope / equals to at z = 1.67. 



Here V denotes for the LDM potentiaL Positive / and negative / correspond to the first 
case and the second case, respectively. 

The reasons why the potential as a function of e shows different features are displayed in 
Fig. [2j The potentials in this figure are referred to the potential at the spherical compound 
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state. At small relative distances z (the red dashed line), reducing e from e = 1 means 
decreasing the surface and thus the corresponding surface potential, while the Coulomb 
potential becomes larger (Fig. |2]^b)). Since the surface potential reduces faster than the 
increasing of the Coulomb potential, the total potential becomes smaller with reducing e, 
as shown by the red dashed line in Fig. [2](a). When the distance between the projectile and 
target increases (black solid line in Fig. |2]) the decreasing speed of the surface area becomes 
smaller. At a critical distance (for instance Zc = 1.67 for ^^°Pd + ^^"Pd), the decreasing of 
the surface potential equals the increasing of the Coulomb potential with reducing e around 
£ = 1, so that the slope at e = 1 becomes zero, as shown by the solid line in Fig. [2]^a)(b). 
While at larger z, shown by the blue dash-dotted line, the surface area is stable with reducing 
e around e = 1 and thus the total potential is increasing while neck starts to evolve from 

€ = 1. 
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FIG. 3: The fusion reactions are classified as hindered (outside the lines) or not-hindered (inside 
the lines) reactions in the neck evolution. On the lines: j de\(^^Y,z=z(i) — 0- The Z and N are 
the proton number and neutron number of the compound nucleus, respectively. See text for more 
details. 

In fact we only need to know the situation at contact [z = zq). By simply checking 
the value of f at z = zq, we know whether the evolution of the neck degree of freedom is 
hindered or not. With this method, different combinations of the Ap and At for different 
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neutron number A^(= A^p + A^t) and proton number Z{= Zp + Zt) of a compound nucleus are 
systematically studied and the results are shown in Fig. [31 The borderlines, corresponding 
to different a, are plotted connecting the reactions for which / equals at e = 1 and z = zq. 
The line of /3-stability is also plotted. Inside the borderlines, the slope / is greater than 
and the reactions are not hindered, while outside the lines it is the opposite. 

Comparing the area encircled with different lines, it also shows that the number of reac- 
tions without neck fusion hindrance for small a is smaller than for large a. This is due to 
the fact that more symmetric reactions correspond to larger Coulomb repulsion force which 
leads the Coulomb potential to increase more quickly with reducing e. Since the encircled 
area is very large even for a = 0, the neck evolution is not hindered for most of the realistic 
reactions. 



IV. FUSION HINDRANCE IN THE RADIAL EVOLUTION 



A. Average of the neck parameter 

We have just concluded that the neck evolution is not hindered for most of the realistic 
reactions. In fact, once at contact a strong potential slope drives the neck degree of freedom 
toward a mono-nucleus shape. On the other hand, the radial and mass-asymmetry degrees 
of freedom face a potential barrier for systems with hindered fusion. As a consequence, the 
neck degree of freedom evolves faster than the other two. It has reached equilibrium before 
the diffusion over the potential barrier has started 1^18|. As most of the di- nucleus evolve 



along z and a direction with an equilibrated neck parameter, we study the fusion hindrance 
in the rest two degrees of freedom at the equilibrated value of e. 

The neck parameter at equilibrium can be determined through the average of e via 

J w[e)de 

where w{e) = e"^*^^-'''^, and T is the temperature of the system. As shown in Fig. HJ {e) 
differs slightly for different reactions. In most cases, (e) is close to 0.1. Therefore we take 
£ = 0.1 in the next calculations. This value is also used in the calculation of fusion cross 
section for ^^Ca induced reactions in the two-step model, and with which the residue cross 
sections are calculated combining the statistical evaporation model [lo|. By introducing 
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FIG. 4: LDM potential as a function of e for different symmetric or asymmetric reactions to form 



a shell energy reduction factor, the calculated residue cross sections show good agreement 
with experimental data from Ca + U to Ca + Bk. 

B. Radial fusion hindrance 

Based on the LDM, when the projectile and target evolve from the contact point to the 
compound state, the parameters {a, z} change from {ao, zq] to {a = 0, z = 0} in the 
same time. However because the LDM potential as a function of a at fixed z and e is 
similar to a parabolic potential centered at a = 0, the evolving system is not hindered in the 
mass-asymmetry direction. Therefore we only need to consider the fusion hindrance in the 
radial direction. Since the friction in the a and z direction is close to each other, the radial 
direction cannot be de-coupled from the mass-asymmetry direction. Thus in the analysis of 
the radial fusion hindrance we will consider a and z direction in the same footing. 

It is easier to understand the radial fusion hindrance in symmetric reactions by simply 
comparing the contact point and the saddle point, shown in Fig. [5l In this figure e has been 
fixed to 0.1. For reactions with light nuclei as ^°Zr + ^''Zr and ^°°Mo + ^'^''Mo, the saddle 
points are located at the right hand side of the contact line {zq= 1.587), which means that 
projectile and target at contact will evolve towards each other automatically by the driving 
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FIG. 5: The relation between the contact point and the saddle point for different reactions. The 
abscissa and ordinate stands for the distance between the two touching nuclei and LDM potential 
(e = 0.1), respectively. The vertical dashed line corresponds to the contact point z = 1.587 for 
symmetric reactions. 



force of the LDM potential and finally form the compound state. While for the reaction 
^^°Pd + ^^°Pd, it is opposite: when the two touching ^^'^Pd reach to the neck equilibrium 
(point "C" in Fig. |5]), they need to overcome an extra-barrier to reach the compound state. 
In other words, the reaction is hindered, or the reaction has radial fusion hindrance. This 
is another fusion hindrance at the contact state besides the neck fusion hindrance. 

However this method cannot be applied to asymmetric reactions because the saddle point 
changes into a ridgeline. But similar to the symmetric case, as we can imagine, if the contact 
point is located inside (left hand side of) the ridge line the reaction will reach the compound 
state without radial fusion hindrance, and v. v. Unfortunately, it is not easy to get the 
ridgeline from a potential surface. To find the fusion path by diffusion, one has to solve 
a Langevin equation as shown in Eq. (6) of Ref. Here, we employ another method: 
switching off the random force, setting the initial momentum in z and a direction as zero 
and then studying the evolution track of the reaction by solving the Newtonian equation. 
The friction tensor is calculated with the wall-and-window formula js]. If the system reaches 
the compound state it means that the contact point is located inside the ridgeline, and thus 
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FIG. 6: Borderlines of the radial fusion hindrance. Reactions located inside the lines do not have 
radial fusion hindrance, while reactions located outside the lines are hindered in the radial direction. 
Three reactions and the line of /3-stability are also plotted. The abscissa and ordinate stands for 
the neutron number N and proton number Z of the compound nucleus, respectively. 

the system is not hindered. Contrarily the contact point is located outside the ridgeline and 
extra-push energy is necessary to form a compound nucleus. In fact this method can also 
be used for the symmetric case. 

With this method we study the radial fusion hindrance in symmetric and asymmetric 
reactions. Fig. [6] shows the location of the borderline between hindered and non-hindered 
reactions for various asymmetries a. For the reactions located on the borderlines, the contact 
points are just located on the ridgeline, while for the reactions located inside the borderlines 
and outside the borderlines, the contact points are located inside and outside of the ridgeline, 
respectively. Therefore all of the reactions located outside of the ridgelines need extra-push 
energy to reach the compound state, or in other words, these reactions have radial fusion 
hindrance. 

Fig. |6] also shows that the number of reactions without fusion hindrance at large a is 
larger than that at small a as it is well known experimentally. This can be easily understood: 
more symmetric reactions (with small a) have larger repulsive Coulomb potential (oc 1 — a^) 
and the ridgeline is closer to or even shifted into the left side of the contact points. In the 
same time, with reducing a, the contact point zq becomes larger, as indicated in Eq. @. 
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This further reduces the number of reactions without fusion hindrance at more symmetric 
reactions. 

It is commonly known from experiment that the reactions ^°Zr + ^°Zr, ^°°Mo + ^"^Mo 
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-|2l|. The 



are not hindered in the fusion reaction while the reaction ^^'^Pd + ^^'^Pd is 
three reactions are also plotted in Fig. [HI Zr+Zr and Mo+Mo are located inside the line of 
a = while Pd+Pd is located outside, showing the agreement with theoretical analysis. It 
is interesting to see that the lower two dots (^^''Hg: Zr+Zr and ^'^^Po: Mo+Mo) is located 
inside all of the lines, showing that all of the reactions to form ^^°Hg and ^°°Po are not 
hindered. While the reactions to form ^sou ^^j^g ^^^^ marked by "^^'^Pd+^^'^Pd") is 
different: because compound nucleus is located outside of the line with a < 0.7 and inside 
of the line with a > 0.7, then the reactions with a < 0.7 are hindered while the reactions 
with a > 0.7 are not. 

Fig. |6]also shows that the number of reactions without fusion hindrance is limited for each 
asymmetric parameter a, and each line for the corresponding a has a .^max- For example, 



Zmax = 87 for a = 0.0 and 2'max = 90 for a = 0.6. According to the extra-push theory [22|. 
■Z^max = 86 for a = and Zmax = 110 for a = 0.6 if Xth = 0.7. The comparison shows that 
the realistic calculations give smaller ^max for larger a, while for the reactions with smaller 
a, the two models give similar Zmax- 

Comparing the Figs. |3] and O it is also found that the encircled area in Fig. |6] is much 
smaller than the area in Fig. [HI and the former one is completely included in the latter one. 
This indicates that the reactions without radial fusion hindrance are apparently without 
neck fusion hindrance, while the reactions with neck fusion hindrance have obviously radial 
fusion hindrance. 

T he sy mmetric case can be compared with the parameterized LDM proposed by Swiatecki 
et al. 22|. In Ref. |22|, the proton number Z and the mass number A of the compound 
nucleus, for which the contact point of the symmetric reaction is located on the saddle point 
of the parameterized LDM, satisfies 

676.48 - 38.534^/172.88 - A/2 ^ ^ ^ 676.48 + 38.534^/172.88 -^/2 

5 + 1352.94/A - - 5 + 1352.94/A ' ^ ' 

Plotting together the above Eq. ([8]) and the borderline at a = of Fig. [6l we get Fig. [71 
It is clear that the lower limit of for fixed Z is almost the same between two models, 
but the upper limit of N is quite different: the parameterized LDM gives much larger iVmax 
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FIG. 7: Comparison of the radial fusion hindrance between the two-center LDM model (solid line) 
and the parameterized LDM model by Swiatecki et al. (dashed line) at a = 0. 



than the two-center LDM model. This is probably caused by the too few parameters used 
to describe the so complicated liquid drop energy in Ref. [22 1. 



V. CONCLUSION 



Fusion reactions between two colliding heavy ions may be hindered or not. This is 
explained by the appearance of an inner barrier after the contact point that hampers the 
formation of a compound nucleus. The liquid drop potential landscape described with the 
two-center parameterization for the deformation allows us to study the frontier between 
hindered and non-hindered reactions as shown in Fig. [3] and El The neck and radial degrees 
of freedom may both have to overcome an inner barrier to reach the compound shape. For 
each of them, we can draw a borderline that classifies reactions with respect to the hindrance. 

It appears that for available nuclei, the neck degree of freedom is not hindered. A strong 
potential slope always drives it toward a mono-nucleus shape. It is not the case for the 
radial degree of freedom. We obtained a borderline that is far more restrictive, i.e. the 
hindrance to the fusion appears for compound nuclei having a charge larger than Z = 87 
when they are formed by a symmetric reaction. The frontier is located at larger values of Z 
for asymmetric reactions, as it is observed experimentally. 

The exact location of the border between hindered and non-hindered reactions might 
depend on the parameterization. It would be therefore very useful to have some experimental 
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results to assess the models. 

Finally, it is important to note that in this paper, we have just studied the appearance 
of the hindrance to the fusion. The amplitude of the hindrance, i.e. the extra-energy that 
is necessar y t o the system to fuse, also depends on dynamical parameters as the friction 
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